A new limit on the Ultra-High-Energy Cosmic-Ray flux with the Westerbork 

Synthesis Radio Telescope 
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A particle cascade (shower) in a dielectric, for example as initiated by an ultra-high energy 
cosmic ray, will have an excess of electrons which will emit coherent Cerenkov radiation, known as 
the Askaryan effect. In this work we study the case in which such a particle shower occurs in a 
medium just below its surface. We show, for the first time, that the radiation transmitted through 
the surface is independent of the depth of the shower below the surface when observed from far 
away, apart from trivial absorption effects. As a direct application we use the recent results of the 
NuMoon project, where a limit on the neutrino flux for energies above 10 22 eV was set using the 
Westerbork Synthesis Radio Telescope by measuring pulsed radio emission from the Moon, to set a 
limit on the flux of ultra-high-energy cosmic rays. 
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I. INTRODUCTION 



Cosmic rays, charged particles traveling through the 
universe, have been observed at energies ranging from 
GeV to above 10 20 eV The highest observed energies 
lie well above energies that can be reached by particle 
accelerators on Earth. These particles are of particular 
interest to astrophysics and particle physics to answer 
fundamental questions about cosmic acceleration mech- 
anisms and particle interactions. The Astrophysics in- 
terest stems from the fact that ultra-high-energy (UHE) 
cosmic rays, particles with energies above 6- 10 19 eV have 
two important properties. They are less deflected by the 
galactic magnetic fields and therefore contain more in- 
formation about their source of origin which is still 
unknown, ft should be noted that the real energy thresh- 
old at which this becomes apparent is dependent on the 
type of particle and scales with their rigidity. In addi- 
tion their sources are not too far from Earth since they 
interact with the cosmic microwave background to pro- 
duce pions and therefore experience substantial energy 
loss over distances of the order of 50 Mpc. This is known 
as the GZK-effect To detect these particles requires 
a large collecting area, because the flux at 60 EeV is only 
l/km 2 /century and drops off faster than E -2 6 . This 
leads to detectors like the Pierre Auger Observatory with 
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a total collecting area of 3000 km 2 Q. However, to de- 
tect particles above 10 21 eV a thousand-fold increase in 
collecting area is expected to be needed. 

As first proposed in Ref. [j| the Moon is a suitable 
candidate with an area of 10 7 km 2 . Detection is based 
on the fact that when a high-energy particle interacts an 
avalanche reaction occurs creating a cascade (shower) of 
particles. The number of particles near the maximum is 
roughly proportional to the total energy in the avalanche 
and is of the order of 10 12 for 10 21 eV @. In this par- 
ticle shower there will be a net excess of electrons due 
to the knock-out of atomic electrons by shower positrons 
and high energy photons (Compton scattering). Simu- 
lations show that this excess amounts to about 20% 0, 
consistent with experimental observations Q. Since all 
particles move with almost the light velocity they are 
closely bunched in the longitudinal as well as the lateral 
direction. In a dielectric this will result in the emission 
of coherent Cerenkov radiation at wavelengths that are 
larger than the typical dimensions of the charge cloud; 
this is known as the Askaryan effect @ . For materials like 
ice, salt or lunar regolith this implies coherent radiation 
at frequencies of 3 GHz and less |7| and this mechanism is 
used in several experiments to detect high-energy cosmic 
neutrinos. Well known examples are the ANITA 
GLUE [HI, LUNASKA [H|, and NuMoon [l]| experi- 
ments. 

For a particle shower deep inside the dielectric, such as 
is usually the case for neutrino-induced showers, the pic- 
ture for the Askaryan effect clearly applies. The subject 
of the present work is to investigate the emitted radiation 
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for the case that the particle shower occurs very close to 
the boundary between the dielectric and vacuum. In this 
case one can imagine that a dielectric layer of minimal 
thickness between the shower and vacuum is required for 
Cerenkov emission to occur. This relates to the concept 
of a formation zone/formation length for Cerenkov radi- 
ation as was introduced in Ref. [14[ and used in calcula- 
tions of the acceptance for detecting high-energy cosmic 
particles [l5| . 

Since neutrinos are weakly interacting particles they 
can traverse many hundreds of kilometers in dense mate- 
rials before interacting. For neutrino-induced showers the 
issue of a possible formation length is thus not essential 
since they typically interact deep in the Moon, compared 
to which any formation length would be negligible. Cos- 
mic rays are particles interacting via the strong interac- 
tion and thus induce a reaction well within a meter from 
the surface of a dense material. For cosmic-ray induced 
showers the possible existence of a formation length for 
the emission of electromagnetic radiation could make the 
difference between being able to observe the shower or 
not especially at wave lengths of the order of a few me- 
ters. The study we present in this paper indicates that 
for the motion of charged particles inside a dielectric the 
concept of a formation zone does not apply. We will 
show that due to the finiteness of a particle trajectory 
inside a dense medium the radiation detected in vacuum 
is independent of the depth of the trajectory below the 
surface, other than for absorption in the medium. This 
situation is very different from that for the opposite ge- 
ometry, an effectively infinite electron beam in vacuum 
inducing Cerenkov emission in a dielectric, as has been 
investigated theoretically in Ref. [l6| and confirmed in 
experiments [l?}- The essence of the difference lies in 
the fact that for particle showers in a medium the track 
length is necessarily finite for which case Cerenkov and 
bremsstrahlung emission cannot really be distinguished, 
as was already noted in Refs. fiH Il9j. 

Using the new finding of an absence of a formation 
length for radiation emitted by a shower in a dense ma- 
terial we have calculated the detection probability for 
observing radio emission from cosmic-ray impacts on the 
lunar surface. The idea to observe this type of emission 
from the Moon with radio telescopes was first proposed 
by [H and the first experimental endeavors in this di- 
rection were carried out with the Parkes telescope [20| . 
It was shown in Ref. [2l[ that observations in the fre- 
quency range of 100-200 MHz (as suggested by Ref. [22|]) 
maximizes the detection probability with small loss of 
sensitivity. This was used in recent lunar observations 
with the Westerbork Synthesis Radio Telescope (WSRT) 
to set a new limit on the flux of ultra-high energy (UHE) 
neutrinos [l3j . Combining these observational data and 
our present calculation of the detection efficiency we set a 
new limit on the flux of cosmic rays at energies in excess 
of 10 22 eV. Such a limit is of interest since recent results 
of the PAO indicate a steepening in the cosmic ray spec- 
trum at the Greisen-Zatscpin-Kuzmin (GZK) energy of 



6 • 10 19 cV [23| . The flux of cosmic rays above this energy 
will be a clear indication of nearby sources and the na- 
ture thereof. In future observations with new generation, 
more sensitive, synthesis telescopes this method may be 
used to measure the flux at GZK energies. 



II. SHALLOW SHOWERS 

The argument that a formation-zone effect will reduce 
the emission for shallow showers hinges on two intuitive 
points: one, that as the distance from the cascade to the 
surface decreases, the emission increasingly becomes in- 
sensitive to the presence of the dielectric; and two, that 
the Cerenkov radiation from such a cascade in vacuum (a 
gedankenexperiment) is zero. In this section, we proceed 
to disprove the second of these two points, and establish 
that current calculational methods, as used in programs 
to estimate the radiated intensity from particle show- 
ers in the Moon, correctly model shallow showers. As 
discussed at the end of this section, this is intuitively 
understandable, because for a finite track the actual ra- 
diation is mostly due to acceleration and deceleration of 
the net charge at the end-points [HI, [l9| . In this process 
coherent radiation is emitted, known as the Askaryan ef- 
fect 1. As is well known [l8|: [l!| for finite trajectories 
one cannot distinguish between Cerenkov radiation and 
Bremsstrahlung. 

We will first (in Section III A[) outline the standard ap- 
proach used in the calculation of the emitted Cerenkov- 
radiation intensity for neutrino induced showers, and in 
Section IIIBI describe how this has been used to calcu- 
late the acceptance for the GLUE [n|, LUNASKA [3, 
NuMoon [l3|, and other experiments. In Section III CI 
we show that this double far-field approach (shower far 
below the surface and the observer far away from the sur- 
face) predicts that a deep cascade (excluding absorption) 
will produce the same observable radiation as a shower 
developing in the vacuum above the lunar surface. Since 
a near-surface cascade (i.e. a cosmic ray induced shower) 
is intermediary to these two extremes one thus expects 
the same radiation as from a deep cascade barring ab- 
sorbtion effects. To show this at a more rigorous level we 
present in Section IIIII an exact treatment which verifies 
that the current double far-field approximation gives the 
correct results even for shallow showers. 



A. Cerenkov radiation from finite particle tracks 

The calculation of the Cerenkov radiation arising from 
finite particle tracks is usually performed assuming the 
observer is in the far-field in an infinite dielectric medium. 
According to classic electromagnetic theory a shower 
with a current density of J z (r,t) = cq(z)S(r — ctz) and 
Jx = Jy = leads in the far field (distance much greater 
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FIG. 1: Geometry of a particle cascade in the Moon. The 
'double-far-field' approach assumes fl>d>L. 



than the shower dimension) to a vector potential 



A z {R,u>) 
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q ^ e -iz U (nco S e a --L)/c dz ^ ^ 



where 9 S is the angle between R and z as shown in 
Fig. (JT|), and n is the index of refraction of the medium. 
The electric field can be calculated to be E = luA z sin 9 S 
where E lies in the plane spanned by the shower and R, 
thus 



E(H, lo) =cj sin 9 



^iujriR/ c 



2-kcR 

q(z)e~ lz ^ ncos ^- 1 ^ c dz . (2) 



The function E(R, u>) can be calculated analytically for 
simple shower profiles @, HH, HH, [25j] , or can be taken from 
parameterisations of detailed simulation results @. 



B. Treatment for deep cascades 

Current methods to calculate the intensity of radio 
emission from Lunar showers use a double-far-field ap- 
proximation. Using the diagram in Fig. ([T]), this means 
that the shower length L is taken to be small compared 
to the shower depth, a, and the travel distance in the 
medium, d, and that the observing distance, R, is much 
larger than the shower depth. The radiation incident 
on the lunar surface can be taken as its far-field (i.e. in 
an infinite uniform medium) solution. This radiation is 
transmitted through the surface according to simple re- 
fraction laws and the strength at the Earth is calculated 
assuming the Earth- Moon distance dominates the 1/R 
term. Thus the observed radiation can be expressed as 



E ohs (R, B a ,u) = t(0i, n) E m (R, 6„,u) e~ d l^ 



(3) 



where t is the transmission coefficient, which depends 
on the incident angle and Lunar refractive index, and 
E m , given by Eq. ©, is the infinite- medium solution of 
the emitted radiation for the Moon's material properties. 
R is taken to be the Earth-Moon distance of approxi- 
mately 3.884x 10 8 m, and the exponential accounts for at- 
tenuation in the medium using the frequency-dependent 
attenuation- length £(uj). 
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FIG. 2: [color online] Radiation pattern as function of ob- 
server angle o at 150 MHz; red drawn curve: a cascade in 
infinite uniform regolith; green dotted: a deep regolith cas- 
cade refracted to the vacuum; black crosses: the total ob- 
served emission from a cascade in vacuum immediately above 
a dielectric boundary. The 'direct' (orange) and 'reflected' 
(blue) components of the vacuum emission are also plotted 
separately. The refracted regolith and total vacuum cases are 
identical. 



The result of Eq. © is plotted in Fig. © func- 
tion of the observation angle 9 Q . The calculation is for 
the plane containing the shower axis and the surface nor- 
mal, so all radiation has a parallel polarisation. Also, we 
take the case of a cascade parallel to the surface (a = 0), 
so Snell's law becomes ncos9 s = cos9 Q . We use Eq. @ 
to calculate E m for a constant charge excess q = — 10 12 e 
moving over a distance L = 3.0 m at velocity v = c. The 
observation frequency is taken as v = 150 MHz, and the 
regolith refractive index as n = 1.8. For waves diverg- 
ing from a point-source, the transmission coefficient for 
parallel polarisation is 



, ( , , 2sin0 o 

n sin 9 a + sin 9 S 

since for a = we have sin 9 S = cos 9i . 

C. Emission in a vacuum 



(4) 



The double-far-field treatment of Section III Bl may be 
expected to break down as the distance d of the cascade 
to the surface becomes small. To place a simple limit on 
near-surface effects, consider the following extreme case: 
a cascade developing immediately above the surface. In 
this case, the radiation seen by an observer will be that 
produced in vacuum (n = 1), but consist of both direct 
and reflected components, with zero path difference due 
to the proximity of the interface. For a = the radiated 
electric field is thus 



E obs (R,e o ,uj) = (l + r{9 0l n))E v (9 ,co) 



(5) 
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where E v is the radiation expected from a cascade in an 
infinite vacuum, and r is the Fresnel reflection coefficient 
(identical for plane and spherical waves). Using Eq. ([2]) 
for E v , radiation of some magnitude will be emitted. In 
Fig. ([2]) the result from Eq. ([5]) is shown using the same 
parameters as used in the double-far-field approximation 
in Section lllBI For the sake of clarity we also plot the di- 
rect and reflected contributions from Eq. ([5]) separately, 
which add to give exactly the same radiation pattern as 
that from Eq. © for d = 0. That is, the emission seen by 
an observer from a cascade immediately above a dielec- 
tric boundary is exactly the same as that from a cascade 
immediately below the surface, which in turn is identical 
to that from a deep cascade if absorption in the medium 
is been ignored. It can be shown analytically in the case 
a = that the two equations give identical results for an 
arbitrary mix of perpendicular and parallel polarisations. 

It seems counter-intuitive that any Cercnkov radiation 
is viewed from a cascade parallel to the surface, since 
radiation emitted at the Cerenkov angle will be totally 
internally reflected. The solution to this apparent contra- 
diction lies in the fact that for a shower of finite length the 
emitted radiation has a large spread around the Cerenkov 
angle (see for example Ref . [HI, ) . For a shower in a 
medium only the radiation emitted at angles larger than 
the Cerenkov angle will penetrate the surface while for 
a theoretical shower in vacuum the Cerenkov angle lies 
at zero degrees. Alternatively one may regard a shower 
of finite extent as corresponding to the acceleration and 
deceleration of charge at the beginning and end of the 
shower [lj| 0j| (equal to the appearance and disappear- 
ance of a moving charge), a picture which has also been 
used to explain the emission of electromagnetic radiation 
of showers induced in air (26j . 

These results indicate that there will be no change 
in the observed radiation as the particle distribution in- 
duced by a UHE particle interaction lies close to the sur- 
face. In the following section this is shown using a more 
rigorous method. 



III. EXACT CALCULATION FOR 
NEAR-SURFACE SHOWERS 

For shallow showers, such as those from UHE cosmic- 
ray interactions with the lunar surface, no far field ap- 
proximation can be made for the radiation reaching the 
surface. For the calculation in the previous section this 
was an essential assumption. We avoid making this ap- 
proximation by performing a complete wave-equation cal- 
culation. 

In the following derivation we consider two half-spaces 
divided by the plane x = 0. For x > the refraction in- 
dex is n' (=1 for vacuum). For x < the refraction index 
is n (=1.8 for the moon). In the lower half plane a parti- 
cle with charge Q and velocity /3 moves from z = —L/2 
to z = L/2 at x = —a, y = 0, passing through z = at 
t = which corresponds to the same geometry as studied 



in the previous section. The four- vector potential for this 
system is determined by the Maxwell equations. Care 
should be taken since the index of refraction is different 
for the two sides of the boundary. 

In the present discussion it is sometimes easier to work 
in the space-time domain, sometimes with energy and 
momentum. The relation between the two is given by 
the usual Fourier transformation 



A(r,t) 



d 3 kdui 
4tt 2 



A(k,w) 



.ik-r— iu)t 



(6) 



In a homogeneous medium with refraction index n the 
vector potential can be calculated [27j from 
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n"to" 47T 
\k l j- * (k,w) = -p M 

c n 

[fc 2 -^]A(k, w ) = ^J(k, w ) 
&■ c 



(7) 



For the problem under consideration the charge density 
is 



=q{z)8(r- vt-a) 



and the current density is 



J(r,i) =vp(r,t) , 



(8) 



(9) 



with v the velocity of the charge, a = —ax the distance 
under the surface, and q(z) the charge of the particle. 
The vector potential can now be written as 



A(k, W ) = 4^ (k ' w) 



cfc 2 



with 



p(k,w) 



q e ik * a sin(±(k z -f c )) 
pc 2^{k z - 



(10) 



(11) 



Note that because v is only non-zero in the z-dircction, 
A z is the only non- vanishing component of A. 

The transmission coefficient for the vector potential 
can be derived from the field equations across the bound- 
ary. In general this is not easy, however, for the special 
case of interest here, this can be done by considering 
the electric field in the y = plane. For electric fields 
generated by a vector potential of the form Ai (r, t) = 
^(Kx+kyv+ktz-ut) wc obtain E, = d a A t - d l A , and 

thus =E x x-\\+E z z-\\=i^^-Al-i^!^-A°. For the 

II x 11 z 11 c k % c k z 

field of the transmitted wave we may write E'n = t| | £7| | . 
In our treatment it is sufficient to calculate A' z since only 
A z is non- vanishing, 



A° 



E ll k ' 



kEu K 



t 



E\\ 



(12) 



where t^n = 2cos£?i/(n'/ncos#i + cos6 t ) (sec Fig. (HJ) 
for the definition of the angles) is the transmission coef- 
ficient of the electric field. At this point wc have made 
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the implicit assumption that the observer is far from the 
surface and that the outgoing waves can be treated as 
plane waves. We still perform a complete integral over 
all waves leading from the source to the surface. For 
transmission parallel to the surface we thus derive that 



2kl 



k x k' x 



rk' 2 



(13) 



where k 2 



ky — k 2 z and k' x = k xX /l 



(n 2 -n' 2 )u 
c 2 k 2 



An incoming wave e lk r - ! "* generates for x > a 
transmitted wave f^Jie r_lw * where ii is given by 
Eq. (|13p . The transmitted radiation (x > 0) can now be 
expressed as 



d 3 k du) 



2kl 



4^ 2 kzK + frk* 



A z (k,uj)e 



ik' -r—iuJt 



(14) 

where A z (k, to) is given by Eq. ([TU|) and k' is a function 
of k. 



To evaluate this expression, first a change of vari- 
ables is made to wave vectors in the n' medium, k' = 



h'2 



72— 



w 2 , k' y = ky, k' z = k z . Subsequently u> is in- 
tegrated over by extending the integral into the complex 
plane where the contributions of the poles of 1 , a 5 

k' 2 — — — 

in A z (k,uj) have to be considered. The poles are cal- 
culated by adding an infinitesimal imaginary part to ac- 
count for causal propagation. This leads to the limitation 
j3ct > L/2. As the next step the d 3 k' integral is written 
as dk' sin 9 d9d<fi where k' x — k' cos 9, k' y = k' sin 6 sin <j> 
and k' z = k' sin 9 cos <f>. The integrals can now be reduced 
to two integrals of the form dk' e lk f^ e ^\ After in- 
tegrating over k' this gives a delta function which can be 
used to perform the integral over <f>. The field can now 
be expressed as 



A z (r,t)=A^(r,t)+A^(r,t) 



(15) 



where the first contribution is 
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W 2 {9) 



2W(8) (g-H 
n'0 



- y 2 sin 2 i 



-T A Re 



rf sin 2 9 -W 2 (9) 



(16) 



and the second contribution is 

Q 



A z 2 \v,t) 



sin 9 d9 - 



W(6)(z -#) + ^ 



W 2 {9) 



2W(8)( z -±) fj 2 2 2 

— wp h y s 



-T 4 Re 



fj 2 sm 2 9-W 2 {9)_ 



(17) 



where 



A. Checking limiting cases 



W{9) = x cos 9 + a cos 6d 1 + 



W{9) = x cos 9 + a cos9\ 1 
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cos 2 9 
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, i L 
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£y cosf 



(18) 



Note that only the real part of the square root contributes 
because of the restriction W(9) < rj sin 9 which is im- 
posed by the (5-function in the derivation. 



In Fig. the vector potential is plotted in the y — 
plane for different times shortly after creation for a ho- 
mogeneous medium. It shows an outgoing pulse traveling 
at the light velocity which is stronger and of shorter du- 
ration in the direction of the Cerenkov cone. At angles 
away from the Cerenkov angle the waves emitted from 
different parts for the current distribution no longer ar- 
rive at the same time. This is reflected in Fig. (|3|) by a 
broadening of the pulse which has the immediate conse- 
quence that the signal is coherent only for lower frequen- 
cies. 

When we include refraction into vacuum the outgoing 
wave has a different structure as shown in Fig. (QJ. Be- 
cause of the refraction at the surface the Cerenkov angle 
is projected at zero degrees, however, there is still radia- 
tion transmitted through the surface. 

The electric- field can be calculated from the vector po- 
tential using E(R,cj) = uj sin 9 (A ( P {K,lu) + a! 2) (R,u;)). 
The Fourier transform is calculated numerically by taking 
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FIG. 3: [color online] The vector potential A z for y = in 
the near field at different times ct=5, 9, 13, 17 m; n = n'=1.8, 
a=0, and L—4 as a function of x and z. The absence of field 
near z = 2 is due to numerical difficulties. 
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FIG. 5: [color online] Electric field strength as a function of 
angle at 150 MHz and 1 GHz at a distance of 1400 km for a 
shower length of 3m, and a charge of Q — 10 12 e, calculated 
according to the double far-field analytic (see Section III A|) 
and the exact numerical (see Section iHTjl methods. 
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FIG. 6: [color online] Electric field strength as a function of 
depth, normalized to that at a = 10~ 2 m, for 3 different angles 
at 150 MHz at a distance of 1400 km for a shower length of 
3 m. 



FIG. 4: [color online] A z for n=1.8, n'=1.0 at ct=20m, L=4, 
and a=l. The absence of field near z = 2 is due to numerical 
difficulties. 

the sum Af (R, w) = Y? t L to e iut Af (R, i) where 

Ag 1 "* (R, t) is calculated numerically at each point. For a 
shower far below the surface (a >> L) viewed from far 
away (R >> a) this should give the same result as the 
analytic result Eq. (J3J. In Fig. (JSJ both the double far 
field analytic result and the exact numerical result are 
compared for two different observing frequencies, show- 
ing that both results are practically indistinguishable. 



B. Formation zone effects 

After having checked the calculation in the far-field 
regime we can use it to calculate the shower at distances 



close to the surface. The result of the calculation is shown 
in Fig. © for 3 different angles. This shows that within 
an accuracy of 1.5% the field is the same for showers at 
depths ranging from 1 cm to 1 km when observed from 
sufficiently far away. The error we attribute to the nu- 
merical calculation and as a second order effect in a/R 
and is largest for the calculation for a = 1 km. The 
important implication of this is that for an observer at 
Earth the observed electric field is independent of depth 
below the lunar surface. There is no shallow shower effect 
which supports the conclusion arrived at in Section [II CI 



IV. NEW LIMIT ON THE UHECR FLUX 

As a first application we will use the present result to 
obtain limits on the flux of UHE cosmic rays using the 
results presented in a recent publication of the 'NuMoon' 
observations of the Moon using the Westcrbork Synthc- 
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sis Radio Telescope (WSRT) [lj. The WSRT consists 
of an array of 14 parabolic antennas of 25 m diameter 
on a 2.7 km East- West line. In the observations we used 
the Low Frequency Front Ends (LFFEs) which cover the 
frequency range 115-180 MHz with full polarization sen- 
sitivity. The Pulsar Machine II (PuMa II) backend HH 
can record a maximum bandwidth of 160 MHz, sampled 
as 8 subbands of 20 MHz each. Only 11 of the 12 WSRT 
dishes with equal spacing were used for this experiment. 
In these observations as part of the NuMoon project, 
the radio spectrum was searched for short, nano-second 
pulses emitted from showers induced in the lunar regolith 
by UHE neutrinos. The data allowed a tightening of the 
bounds on the neutrino flux at high energies (T3 |. 

When an UHE neutrino interacts, most of the energy 
is carried away by the emerging lepton which, in general, 
does not produce a detectable signal while only about 
20% of the energy is deposited in a hadronic shower which 
emits a signal that can be detected at Earth. When a cos- 
mic ray impinges on the lunar surface all its energy will 
be converted into a hadronic cascade of energetic par- 
ticles. This cascade will commence right at the surface 
and the shower maximum is thus within meters from the 
surface. Due to the absence of a formation zone these 
events should thus also give a signal of the characteris- 
tics of the ones that were searched for in the observations 
of Ref. jl3| . With its surface area of the order of 10 7 km 2 , 
our finding of the absence of a formation zone thus shows 
that the Moon can be used as a sensitive cosmic-ray de- 
tector. 

As the first step in the WSRT observations the narrow 
band Radio Frequency Interference (RFI) is filtered from 
the recorded time series data for each subband (with a 
sampling frequency of 40 MHz) and the dispersion due 
to the ionosphere of the Earth is corrected for. Short, 
nano-second, pulses emitted from the Moon correspond 
to strong pulses with large bandwidth. To search for 
these pulses 5-time-sample-summed power spectra (so 
called iVspectra) were constructed for all subbands. The 
data were kept for further processing when in all four sub- 
bands beamed at the same side of the Moon a value for 
P5 larger than a certain threshold was found within a 
certain maximum time offset. In a subsequent analysis 
additional constraints were imposed, such as eliminating 
pulses that had a long time duration and pulses that were 
found in both beams within a certain time limit. The 
largest remaining pulse had a strength of S = 152 kJy 
which is a factor 3 larger than expected for pure statis- 
tical noise. To account for dead-time issues and filtering 
inefficiencies a complete simulation of the data taking 
was performed where short pulses with a random time 
offset have been added to the raw data (considered as 
background). The detection efficiency (DE) was deter- 
mined as the fraction of inserted pulses that is retrieved 
after applying all trigger conditions and cuts that were 
used in the analysis. The System Equivalent Flux Den- 
sity (SEFD) for the WSRT, averaged over the frequency 
range under consideration, is a 2 = 400 Jy per time sam- 
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FIG. 7: [color online] The currently established cosmic-ray 
flux limit from WSRT observations [13|] (thick red drawn line) 
is compared to the flux determined by the Pierre Auger Ob- 
servatory [32| (data points with error bars) and a simple poly- 
nomial expansion (black line, see text). Also the prospective 
flux sensitivities are indicated that can be obtained with LO- 
FAR H and SKA H observations. 



pie. 

A more detailed des crip tion of the followed procedure 
can be found in Ref. [13[ where it is concluded that in 
46.7 hours of observation no pulses from the Moon were 
detected with a strength exceeding 240 kJy with a 87.5% 
probability. Simulation calculations have been performed 
to convert this to a flux limit. In the simulation cos- 
mic rays of a certain energy E s hit the lunar surface at 
arbitrary angles and create a particle cascades. Based 
on Monte Carlo simulations the intensity of radio-waves 
emitted from such a cascade in the lunar regolith have 
been parameterized as function of emission angle and fre- 
quency [29M3 1( | . The emitted radiation is passed through 
the lunar surface following the usual laws of wave refrac- 
tion which, due to the absence of formation zone effects, 
also applies to the case where the particle cascade occurs 
just below the surface. From the surface of the Moon 
to Earth the intensity follows the usual inverse square 
law. The details of the simulations are described in de- 
tail in Ref. [2l| . U sing the model- independent procedure 



described in Ref. [24| the non-observation of pulses of 



a certain strength can thus be converted in an energy- 
dependent 90% confidence limit on the cosmic-ray flux 
as shown in Fig. ((7|) by running the simulation for differ- 
ent values of . In the simulations the effects of surface 



roughness can be ignored [21 [ . 

The present limit for the flux of cosmic rays that fol- 
lows from the existing WSRT observations |13( is well 
above what could be expected based on the observa- 
tions made by the Pierre Auger Observatory (32|. The 
thin dotted straight line in Fig. ([7]) shows the model- 
independent differential flux limit, comparable to the 
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presently set limits, for observations of the Pierre Auger 
Observatory, showing that the limit set by the WSRT 
observations is considerably lower albeit at considerably 
higher energies. The thick line corresponds to a p olyno- 
mial expansion E~ 4 - 3 as advocated in Ref. [32[ where 
the grey band corresponds to an uncertainty in the ex- 
ponent of 0.8. The lower limit of the exponent has been 
taken from Ref. (33|. The quoted values for the flux in 
Ref. [33[ lie well above the ones shown in Fig. ([7]) which 
might be due to an uncertainty in the energy calibra- 
tion [32l ]. Future observations with new-generation ra- 
dio telescopes such as LOFAR [34[ or SKA 35 1 should 
reach much higher sensitivity for pulse detection, result- 
ing in correspondingly lower energy thresholds as shown 
in Fig. ([7]). We show the sensitivity that can be reached 
in a one week measurement using the LOFAR telescope 
where the drawn curve uses only the core stations (SEFD 
of 93 Jy, using 50 MHz bandwidth) while the long-dashed 
curve uses all E-LOFAR stations (SEFD of 30 Jy, using 
full bandwidth). We have assumed here a 100% moon 
coverage and a detection threshold of 6cr where a is the 
amplitude of the noise. In a one week measurement with 
the future SKA telescope (SEFD of 1.8 Jy) the results 
depend on the frequencies used for the observation. At 
lower frequencies (100-300 MHz band, SKA-1 in Fig. ©) 
one is sensitive to a smaller flux while at intermediate 
frequencies (300-500 MHz band, SKA-m in Fig. ([7])) one 
is sensitive to cosmic-rays of lower energy. The increased 
sensitivity will make this method sensitive to cosmic ray 
energies of the order of 10 20 eV where, due to the large 
collecting area, competitive measurements are possible. 



V. CONCLUSIONS 

We have shown that the concept of a formation zone 
does not apply for the emission of electromagnetic ra- 



diation from a moving charge distribution over a finite 
distance inside a dielectric emitting Cerenkov radiation. 
In particular we have considered the system where the 
charges move at close proximity to the surface separat- 
ing the dielectric and vacuum. We have shown that the 
radiation penetrating the surface is independent of the 
distance of the charge distribution to the surface even 
for distances that are much smaller than the wavelength, 
provided the observer is sufficiently far away from the 
surface. In principle we have shown the absence of a for- 
mation zone only for a charge distribution with a block 
profile, however, the superposition principle can be used 
to show that this conclusion extents to showers with a 
realistic profile. 

One field where this finding has a large impact is in 
the calculation of the acceptance of large scale cosmic- 
ray and neutrino detectors. As one application we have 
calculated the acceptance of the observations for the Nu- 
Moon project to cosmic rays and used this to derive a 
limit on the flux of cosmic rays for energies in excess of 
2 x 10 22 eV. If there would have been a formation length 
of the order of the wavelength of the observed radiation, 
the acceptance would be vanishingly small. We instead 
find that this approach will offer a competitive means of 
detecting the flux of cosmic rays with energies in excess 
of 10 20 eV with future synthesis radio telescopes. 
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